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Abstract 

Novel chirally symmetric fermion actions containing the minimum amount of fermion doubling have been 
recently proposed in the literature. We study the symmetries and renormalization of these actions and find 
that in each case, discrete symmetries, such as parity and time-reversal, are explicitly broken. Consequently, 
when the gauge interactions are included, these theories radiatively generate relevant and marginal operators. 
The restoration of these symmetries and the approach to the continuum limit thus require the fine-tuning of 
several parameters. With some assumptions, we show that this behavior is expected for actions displaying 
minimal fermion doubling. 
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I. INTRODUCTION 



A notorious hurdle in the study of lattice field theories is the phenomenon of fermion doubling. 
In fact, most actions containing only one fermion field per lattice site while retaining at least a 
remnant chiral symmetry describe several fermion species in the continuum limit. The presence of 
these extra degrees of freedom is, to a certain degree, unavoidable. The Nielsen-Ninomiya [H, Q, 0] 
theorem guarantees that under certain very general assumptions, chirally symmetric actions lead to 
a doubling of the number of fermions one would naively expect. Fermion actions which bypass the 
Nielsen-Ninomiya no-go theorem are known but are computationally very expensive. For actions 
possessing hypercubic symmetry there is one fermion doubling for every space time direction and 
thus, in four dimensions, one Dirac fermion on each lattice site leads to sixteen Dirac fermions 
in the continuum limit. Their use in QCD where only two (or three) light fermions exist is 
problematic. In the past, Karsten 0] and Wilczek 0] have suggested four-dimensional actions with 
the minimal number of doublings allowed by the Nielsen-Ninomiya theorem, namely, one doubling, 
corresponding to two fermions in the continuum. Unfortunately, these actions require the fine 
tuning of some parameters to achieve the continuum limit, which is a daunting numerical task, and 
have not been extensively used. The subject was recently revived with the suggestion by Creutz 0] 
of a graphene inspired class of actions, some of them defined on non-hypercubic lattices, possessing 
both chiral symmetry and the minimal fermion doubling. Soon afterwards Borigi 0] suggested a 
related action and clarified some points of Creutz's proposal. 

In this letter, we point out that these actions [1, H, |g, 0] break discrete symmetries like parity 
and time reversal and therefore when gauge interactions are included, one should expect to fine 
tune relevant and marginal operators in order to have a continuum limit. The need for fine tuning, 
of course, hinders the usability of these actions in practical calculations. Before presenting a 
systematic analysis, it is useful to make our main point using the Borigi action as an example, 
which, in momentum space, can be written in the form 



aD B {p) = -^—^ + ^74 ^2 cosiapn) + i 

M=l 3=1 



Sj(ap), (1) 



where 



s 2 (ap) 
K s 3 (ap) / 




fsm(api)\ 

sin(ap 2 ) 
sin(aj> 3 ) 
\sm(ap 4 )J 



(2) 



This action breaks the parity symmetry, under which ip{p,Pi) — > 74'0(— p, p±), charge conjugation, 
'tpiPiPi) ~ C'tjp- (p,Pi), and time-reversal, ip{p,Pi) — ► 7574 V'CP) ~Pi)- Consequently, in the effec- 
tive (Symanzik) action valid at distances larger than the lattice spacing a but smaller than the 
confinement radius, the time-reversal breaking operator 

- tp «74 i> (3) 
a 

appears in the theory. This relevant operator is not part of the continuum limit of QCD and 
therefore requires its explicit inclusion in the action with a finely tuned coefficient. This is an 
example of the fine tuning problem alluded to above. The general analysis of the symmetries and 
the operator content of the Symanzik action of these theories requires more care and is the central 
subject of this paper. We will see that these novel actions require multiple fine tunings of dimension 
three and four operators. 
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We begin with a general analysis of the Borigi and Creutz actions in Sec. HH rewriting these 
actions to expose their flavor structure. This allows us to identify symmetries of the action and 
construct the first few terms of the Symanzik effective continuum theory, which is done in Sec. IIII1 
In Sec. \TV\ we comment on the analysis for the Wilczek action. We conclude in Sec. |V] by arguing 
that any chirally symmetric action of this type with minimal doubling will also require the fine 
tuning of relevant operators. 



II. THE BORigi AND CREUTZ LATTICE ACTIONS 

Here we consider the fermion actions proposed by Borigi [7] and Creutz 0]. These actions can 
be written in the general form 



Sbc 



x,ii 



ip(x) ( Ef, + ) ip(x + n) - ip{x) ( ~ M - iB<j 4 )ip(x - fi) - 2BCip(x)i'j4 i>(x) 



where the matrices are linear combinations of the Dirac matrices 7^ , 

3 

3=1 



(4) 



(5) 



with 



A 



(+1 +1 +1\ 

+1 -1 -1 

-1 -1 +1 

V-l +1 -1/ 



(6) 



The notation for S M is only a shorthand because H M does not transform as a four-vector under the 
spinor representation of the Lorentz group. The Borici-Creutz action depends on two parameters 
B, and C. When B = 1, and C = l/v2, Eq. (|4|) reduces to the Borigi action. 

This action, Eq. (j3j), can be diagonalized by expanding the field ip{x) in Fourier modes. With 



* <£p 



ipip) , 



we find 



Sbc 



7T _ 



(7) 



(8) 



where the propagator Dbc(p) is given by 



Dbc(p) 



Uz-n sinp^ + iBj^cospfj, - C) 



(9) 



From the last term in the propagator it is clear that discrete symmetries C and T are broken. 
Furthermore H4 is odd under parity, while sinp4 is even leading to broken P in Eq. (|9|). 

Provided B 7^ 0, and < C < 1, there are exactly two poles of the propagator corresponding 
to massless fermions, namely at PfP = (p,p,p,p) and py — '"' 



Pfj, , where cosp = C. From the 
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momentum space action, it is straightforward to define flavored quark fields. Momentum modes 
near pjP are described by the field Q^\p), while those near pffi are described by Q^ 2 \p)- The 
exact partition of the Brillouin zone into these two non-overlapping regions is irrelevant. Q^ip) 
should have support in a momentum hypersphere centered at and Q^ 2 \p) should have support 
in a hypersphere about p^ . 1 The union of the two partitions must equal the full Brillouin zone. 
In this scheme, the action can be decomposed into two terms 



Sbc 



/7T 
-7T 



* d*p 



Q {1 \p)D-\p)Q { - l \p) + Q [2 \p)D-\p)Q { - 2 \p) 



(10) 



(2vr) 4 

It is convenient to shift both terms so that the momenta are measured from zero, i.e. 

Sbc = f ^[& 1 \<l+tf ) )D-\ q + P ^ 

(11) 

where periodicity of the momentum fields allows the integration to remain over a single Brillouin 
zone. The propagator near has the form 



ah 



D(q±pW) 



(12) 



Above S is given by S = Vl-C 2 . 

We can compactly write the action to expose flavor symmetries with the definition of an isodou- 
blet Q, 



VW-\TQW(q-pW))' 



(13) 



where T is a matrix evoking the similarity transformation: {7^- — > 7^,74 — > —74}- Thus we find 



Sbc 



(2tt)< 



J2 \ Q(q) [(CE^ - BS lA ) ® l] i sin ^ Q(q) 

+ Q(q) UsEp + SC74) ® r 3 l i(cos g M - l)Q(g) 1 . 



(14) 



The isospin matrix r 3 has the familiar definition, r 3 = diag(l, —1). 

At this point, we identify a new set of Dirac matrices 7' given by 7^ = X^7^ a ^> with 



1 

a = - 
2 



(+1 +1 -1 -1\ 
+1 -1 -1 +1 
+ 1 -1 +1 -1 

V-i -1 -1 -1/ 



(15) 



1 Notice that a free valence quark of the first variety, for example, can thus be constructed from an operator of the 
form 

where V is a momentum space hypersphere centered about p^ 1 '. When the theory is gauged, a Wilson line is 
needed for gauge invariance. Calculation of correlation functions of definite flavor requires computations at the 
level of all-to-all propagators. 
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Due to the orthogonality of a^ u , the 7^ is indeed a set of Dirac matrices [7:]. We then rewrite the 
Borici-Creutz action in the form 



Sbc 



-n (2VT) 4 



£ 1 0(g) [ \2C^ --(C-BS)^2 llj ® l] i sm % 0(g) 

+ Q{q) \ ( 25 7 ; - ~(S + 5C) l'u ] ® t 3 1 *(cos g M - l)O(g) 1 • (16) 



One can produce orthogonal axes of the continuum action with the choice C = BS. For B = 3, this 
corresponds to the value C = 3/\/l0 as shown in [6J. For B = 1, this corresponds to C = 1/V2, 
and results in the Borigi action [7j]. For an arbitrary value of B, the choice C = B/y/l + B 2 
produces orthogonal axes. By defining the 7^ matrices, we have absorbed the momentum rotation 
into the Dirac algebra. After rescaling the fields, the orthogonal axis Borici-Creutz action has the 
form 



oorthog 
^BC 



n d 4 q 



^Q(9)(->£®l)» si* g M 0(g) 



+ 1 0(g) 



lu 



-(l + B 2 )Y,l'v)®T 3 ]i(cosq,-l)Q( q ) . (17) 



The first term in Eq. (|17p is hypercubic invariant but the second is not. Thus no additional 
symmetry emerges upon choosing orthogonal axes. 

When C 7^ BS, the lattice hyperplanes do not constitute an orthogonal coordinate system for 
the particle momenta. In the continuum limit, the action Eq. (|16[) tends to ij^k^, where the k^ 
specify the particle momenta in an orthogonal system. This yields ku = R^qu, where as a matrix, 
R^ 1 is given by 



1 

AC 



'+1 +1 +1 
+1 -1 -1 



-1 -1 +1 

v- 1 + 1 - 1 -4sJ 



(18) 



In terms of physical particle momenta, we have 



Sbc 



d 4 k 
(2tt)< 



Q(k) \ 2C% - -(C - BS) J2 ll ® ll i sHR-'k), Q(k) 



+ Q(k) (25 7 ;- \{S + BC)Y,l'A iUo S {R~ l k)^-l)Q{k)\.{\9) 



The symmetries of the properly gauged version of Eq. (|19[) are specifically as follows: 

1. Gauge invariance 

2. Discrete spacetime translation invariance 

3. U(1)b baryon number 

4. U(1)l <8> U(1)r chiral symmetry corresponding to t 3 

5. Combined CPT invariance (but none individually) 
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6. Combined IC invariance, where I is the isospin rotation I = exp (±«7rr 1 /2) 

7. S4 invariance under permutations involving all four hyperplane axes. 

The last symmetry can be translated into fairly complicated transformations among the physical 
axes of particle motion. We have, however, found it easier to work with the original hyperplane 
axes. While there is an axial U(1)a symmetry possessed by Eq. (|16p . near the continuum it will be 
spoiled by the anomaly. Lastly spontaneous symmetry breaking reduces the chiral symmetry down 
to U(l)v flavor symmetry generated by r 3 . It might be possible that additional non-standard 
symmetries emerge for particular values of B and C, such as for the values B = y/b cot(7r/5) and 
C = cos(-7r/5) suggested by analogy with graphene jg]. 



III. SYMANZIK EFFECTIVE THEORY 



To investigate discretization errors, we consider the continuum limit. Using the seven symme- 
tries of the Borigi-Creutz action, it is straightforward to construct the Symanzik effective theory, 
valid at distances larger than the lattice spacing but smaller than the confinement radius 0, 0] . 
The general form of the Symanzik theory is 

00 

71=3 j 

where the sum on n is over the operator dimension and the sum on j is over all operators of a given 
dimension. As the symmetries of the Borigi-Creutz action are non-standard, however, we shall 
begin first with a naive demonstration to produce a term generated from radiative corrections. 

Consider the lattice spacing and gauge coupling sufficiently small. Accordingly we can treat the 
theory as weakly coupled. Moreover, we can expand the action in Eq. (|19p in powers of the lattice 
spacing and perform a Fourier transform to position space. Up to 0(a), the theory becomes (in 
dimensionful units) 

C eS = Q( 7j u ® l)D^Q - -a-Q \2S^ - \{S + BC) £ >A r 3 '(iT 1 D) ^R' 1 D) M Q , 

(21) 

with the gauge covariant derivative given by = 8^ — igA^, where is valued in the algebra of 
the gauge group. The gauge field strength tensor is given through the relation [D^, D„] = —igF^ v . 




FIG. 1: Self-energy tadpole which generates contributions to a relevant operator. 

Now we use the effective action in Eq. (|2ip to investigate the fermion self-energy. Working 
perturbatively in the coupling constant, there are discretization corrections which arise from an 
insertion of the 0(a) dimension five operator C5 = Q(iC f _ l ®T 3 )(R~ 1 D)^(R~ 1 D)^Q. Here £/j is just 
a replacement for the combination of 7^ matrices appearing in Eq. (|2ip . The C M do not in general 
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satisfy the Dirac algebra. A self-energy diagram generated from inserting this operator is shown 
in Figure [TJ This tadpole is straightforward to evaluate. Schematically we find 

^^aa s (iC M ®r 3 )u(p) ~ — u(p) ® t 3 )u( P ), (22) 

q » a i* 

which corresponds to a relevant operator in the Symanzik action of the form 

03 = ]TQ«^t 3 )Q. (23) 

A* 

Now it could be possible that additional 0(a s ) diagrams happen to cancel this 0(a~ 1 ) term. 2 
However, without a symmetry protecting this term from being generated by the interactions, it 
likely receives contributions from all orders in the gauge coupling and appears in the Symanzik 
action with a non-perturbatively determined coefficient. One can verify this operator, O3, satisfies 
the seven symmetries of the Borigi-Creutz lattice action. 

(7) 

Using the seven symmetries of the Borici-Creutz action, we construct all relevant, 0\ , and 

a) 

marginal, 0\ , operators present in the Symanzik effective theory C e Q, 

^= l -H W } + E °V°V + 0(a). (24) 
3 3 

The relevant operators are 

of = g(i 74 ® t 3 )q 

Of = Q( 7 475 ® t 3 )Q. 

One notes that Y2 u 'y' I/ = —274, which makes obvious the permutation symmetry of the original 
hyperplane axes. The operator O3 can be eliminated by a field redefinition, see Sect. IIVI for a 
discussion of the analogous operator for the Wilczek action. 
The marginal operators in the Symanzik theory are 



of 


= Q( 7m (8 l)DnQ 


of 


= Q{li ® 1)D 4 Q 


of 


= Q(«7m75 ® tyDpQ 


of 


= Q(«7475 ® 1)D*Q 


of 




of 




of 


= FinFi/x 


or 


= F^F^, 



2 There are also operators of higher dimension that contribute to the fermion self-energy at 0(a~ 1 ). The reason 
being that higher dimensional operators can produce terms ~ {aq) n , where q is a loop momentum that can 
reach the cutoff. Thus all such terms are order unity. One must use lattice perturbation theory to determine all 
contributions. Our example is meant to demonstrate the existence of such terms, not to determine their coefficients 
rigorously. 
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FIG. 2: High momentum gluon interaction which leads to flavor exchange. Here the gluon momentum is 



2<i 



which is specific to the Borii action. In general, the flavor changing momentum is - cos C. 



where F is the dual tensor to F, namely F^ u = ^nvpa-Fpa- Beyond this order there is an explosion 
in the number of operators. For example, at 0{a) there are irrelevant operators that lead to 
chromoelectric dipole moments. At 0(a 2 ), one has various four quark operators. Analogous 
to the case of taste breaking in staggered fermions, there are flavor changing interactions. We 
exemplify this in Figure where a high momentum gluon exchange leads to mixed-flavor four- 
quark interactions. Notice that the individual flavor numbers are conserved in accordance with the 
U(1)b <8) U(l)v symmetry of the action. 



IV. THE WILCZEK ACTION 



Karsten 0] and Wilczek [j| pointed out, over twenty years ago, the possibility of minimal 
doubling with an explicit construction. The Wilczek action has the form 



2 

X,fJt 



A [i'i^n^ix + j) + i>(x)i~f 4 ii(x - j)] , (25) 



2 

x,3 



with A > i as a parameter. The action is diagonal in momentum space 



2 

"* d A k -~ 



Sw= / ^^(^"^M*)' ( 26 ) 



with the propagator D{k) given by 



D(k) 



i'y^ sin — 174A ^^(cos kj — 1) 



(27) 



There are exactly two poles of the propagator, namely at the value s fc^ = (0,0,0,0), and kffl 



(0,0, 0,7r). As this action has been analyzed in detail previously [lOJ, we shall not pursue the 
decomposition of the fields in terms of flavor and chirality. Instead we will focus on the discrete 
symmetries of the action: 

1. Cubic invariance 

2. CT invariance 

3. P invariance 
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FIG. 3: Two dimensional slice of momentum space. Solid squares denote poles, or potential poles of the 
propagator. 



4. invar iance, 

where is the time link-reflection symmetry operator. For the Borigi and Creutz actions, is not 
a symmetry. 

Given these discrete symmetries, it is straightforward to write down the relevant and marginal 
terms in the Symanzik action. There is one relevant operator 



e> 3 = 



and several marginal operators 



O 
O 
O 
O 



(i) 

4 

(2) 
4 

(3) 
4 

('!) 



As with anisotropic lattices, the breaking of hypercubic symmetry down to cubic symmetry requires 
a speed of light tuning. It has been argued in Ref. 10] that the relevant operator can be absorbed 
by a field redefinition, namely 



exp 



2C3X4 



a(l + c 



(3)> 



(28) 



where C3 is the non-perturbatively determined coefficient of the operator O3 in the Symanzik action, 

(3) (3) 

while C4 that of 0\ . The field redefinition, however, modifies the boundary condition and just 
moves the fine tuning problem from the bulk of the action to the boundary. The physical reality of 
this term can be ascertained by noticing that it plays the role of an imaginary chemical potential. 
The relevant operator thus cannot be avoided in the case of the Wilczek action. 



V. CONCLUSION 

We have investigated above different lattice actions with minimal fermion doubling. In partic- 
ular, we demonstrated that the Borigi-Creutz action has an exact U(1)l <8> U(1)r chiral symmetry. 
For each action, however, there are relevant and marginal operators appearing in the Symanzik 
effective theory. Breaking of discrete symmetries, such as parity and time-reversal, is the culprit 
for the appearance of such operators. The Wilczek action has the least number of such operators. 

A natural question arising from our analysis is whether a chirally symmetric action with minimal 
fermion doubling which does not generate dimension three operators is possible. This would indeed 
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occur if such a putative action preserved PT symmetry as that symmetry is sufficient to forbid the 
dimension three operators listed. It is, however, easy to see that such an action does not exist. To 
have minimal fermion doubling, it is clear that hypercubic invariance must be broken. Heuristically 
one can argue that parity or time-reversal must also be broken from Figure O Here we plot a two- 
dimensional slice of momentum space. A shift in momentum space can always be made so that one 
pole of the propagator is at a corner of the slice depicted. Poles of the propagator are shown in the 
case of ordinary fermion doubling, and two minimal doubling scenarios. Notice that in the absence 
of symmetric doublers, there is always an asymmetry about one axis. Hence, depending on which 
axis is asymmetric, either P, or T is broken. The dashed line shows a hypothetical scenario in 
which the product PT is preserved. If we write the general chirally symmetric momentum space 
action as 



then we require that each f ^ and be odd, e.g.. fn(—p v ) = —f^(Pu), so that PT is preserved. 
Furthermore the action must be local which requires the continuity of /«, g^ and periodicity. The 
only odd functions which are continuous and periodic within the Brillouin zone vanish at the 
boundaries ir/a and —ir/a (as well as at the origin). Thus there is no PT symmetric action which 
has only minimal doubling. We have reasoned that chirally symmetric minimal fermion doubling 
requires broken P, or T, and broken PT. Without these symmetries intact, however, dimension 
three operators are allowed, e.g. ipi'jjip for j spatial, or ^7475 ^ in the case of broken P, and 
-i/>7j75?/> for j spatial, or ipi^ip for the case of broken T. Thus chirally symmetric minimal doubling 
actions, of the form considered here, require the fine tuning of relevant operators. There may 
be additional non-standard symmetries, however, that emerge for particular values of parameters. 
Such symmetries are very interesting because they could potentially eliminate relevant operators. 
It is thus worthwhile to expose such symmetries. 
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